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Abstract 

We study non-Gaussianities in the primordial perturbations in single field inflation where there 
is radiation era prior to inflation. Inflation takes place when the energy density of radiation drops 
below the value of the potential of a coherent scalar fleld. We compute the thermal average of the 
two, three and four point correlation functions of inflaton fluctuations. The three point function 
is proportional to the slow roll parameters and there is an amplification in f^L by a factor of 65 
to 90 due to the contribution of the thermal bath, and we conclude that the bispectrum is in the 
range of detectability with the 21-cm anisotropy measurements. The four point function on the 
other hand appears in this case due to the thermal averaging and the fact that thermal averaging 
of four-point correlation is not the same as the square of the thermal averaging of the two-point 
function. Due to this fact ttvl is not proportional to the slow-roll parameters and can be as large 
as —42. The non-Gaussianities in the four point correlation of the order 10 can also be detected 
by 21-cm background observations. We conclude that a signature of thermal inflatons is a large 
trispectrum non-Gaussianity compared to the bispectrum non-Gaussianity. 
PACS : 98.80.Jk, 98.80.Cq, ll.lO.Wx 
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1. INTRODUCTION 



The experimental determination of signature of primordial non-Gaussianity in the CMB 
spectrum is of great interest as non-Gaussianity can give us insights into the dynamics 
of inflation models Interacting fields have non-Gaussian correlations proportional 

to the coupling strength, which for the infiaton is unmeasurably small j^. In quasi de- 
Sitter space, single field inflation models predict a level of non-Gaussianity proportional to 
the slow roll parameters {4, s], Q]. The current constraint on CMB bispectrum from the 
WMAP Syr data is —151 < f^j^ < 253 (95% CL) j^, /atl being the bispectrum non-linear 
parameter and a measure of non-Gaussianity. The single field inflation models prediction 
of Jnl ~ 10~^ is too small to be detectable in WMAP or the upcoming PLANCK mission 
where non-Gaussianities at the level of Jnl ~ 5 can be probed 8j. The most sensitive 
probe of primordial non-Gaussianities can come from the measurement of anisotropics of 
the 21-cm background whose bispectrum can probe fNL< 0-1 [sl, [l^. WMAP constrains 
the non-Gaussianity from trispectrum at |t7vl| < 10^ [ll| while PLANCK is expected to 
reach the sensitivity upto \tnl\ ~ 560 [l^. The anisotropics of the 21-cm background can 



constrain the trispectrum of primordial perturbations to the level of ~ 10 [lO], which is still 
too large compared to the predictions of the single field models of tnl — {^Inl)"^ [isl, \l4 \. 

It was shown by Gangui et al. [15] and more recently in [g], IigI that if the initial state 
of the inflatons is not the Bunch-Davies vacuum but some excited state then there is an 
enhancement of the non-Gaussianity from such initial state effects. A natural example of a 
non-Bunch Davies initial state arises if there is a pre-inflation radiation era prior to inflation 
Inflation takes place when the energy density of radiation drops below the value of 
the potential of a coherent scalar field. In such models it is seen that the power-spectrum 
is enhanced at low k which can be used to put constraint on the comoving temperature 
at the time of inflation |17I]. Inflation scenario with a pre-radiation era have an interesting 
prediction that the B-mode polarization spectrum is enhanced at low / due the contribution 



of thermal gravitons 18|, Il9 |. 



The scenario of thermal initial condition is very general and would be applicable for any 
model of inflation if there was a pre-inflationary radiation dominated era. The effects of 
the initial thermal era to be observable either in the CMB anisotropy spectrum or in the 
non-gaussianities the perturbations entering the horizon today should have left the de-Sitter 
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horizon at a temperature T not too small compared to H (the Hubble parameter at the time 
of inflation). If there were a large number of e- foldings prior to the present perturbation 
modes leaving the inflation horizon then the effect of the pre-inflationary thermal era would 
be unobservable. In models where the total number of e-foldings are just enough to solve 
the flatness and horizon problems, there can be a imprint of the spatial curvature at the 
time of inflation on the power spectrum. This has been studied in [20]. A natural model 
where inflation commences just as the temperature falls below a critical temperature and 



is of limited duration is where a fermion pair 



inflaton. Such models have been studied in 
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"orms a scalar condensate which acts as the 



22|. 



In this paper we study non-Gaussianities in the primordial perturbation in single field 
inflation where there is radiation era prior to inflation. The thermal background of inflaton, 
gravitons and other fields is decoupled from the actual dynamical evolution of the inflaton 
unlike in the warm inflation models {23], where there can be large non-Gaussianities 24, 25| 
due to dissipative coupling between the inflaton and the radiation bath. In this model 
the temperature of the decoupled radiation bath goes down as Tph = T/a where T is the 
constant comoving temperature. The thermal distribution functions which depend on the 
ratio = ^ (where k is the comoving wavenumber of the perturbations) retain the same 
form during inflation. 

We calculate the thermal average of the three point correlation function, otherwise known 
as the bispectrum, of the comoving density perturbations. The non-linear parameter f^L 
being a measure of non-Gaussianity due to bispectrum is turned out to be a function of 
the magnitude of the three momenta of comoving perturbations. In order to quantify J'nl 
three distinct configurations of those momenta are analyzed namely the "squeezed" triangle, 
the "equilateral" triangle and the "folded" triangle and it is observed that the maximum 
contribution for /jvl comes from the "equilateral" configuration. We find that the thermal 
contribution can result in the enhancement of J'nl by factors ranging from 65 — 90. 

We show that due to the presence of the initial temperature the contribution to four 
point correlation function of the density perturbations comes from different factors other 
than the mere disconnected diagrams. We evaluate the thermal average of the four point 
function and calculate the contribution of the thermal initial states to tj^h. We find that in 
the leading order t^vl is independent of the slow roll parameters and we find that it can be 
as large as tat^ ~ —42. 
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We organize the rest of the paper as follows. In Section [21 we compute the thermal average 
of the two point correlation of the density perturbations. We see that the two point function 
is enhanced at low k due to the contribution of the thermal inflatons as observed in 'y\ . 
In Section [31 we compute the thermal averaged three-point correlations and calculate /tvl 
for various momentum configurations. In Section [H we compute the thermal contribution to 
tnl- We discuss the feasibility of measuring the non-Gaussianity due to thermal inflatons in 
the Concluding section. We have given a outline of the calculation of the non-Gaussianity 
parameter /tvl at zero temperatures for the single field slow- roll inflation models in the 
Appendix. 



2. THERMAL AVERAGE OF INFLATON POWER SPECTRA 

If there was a radiation era prior to inflation one expects a thermal distribution of infla- 
tons to be present which might have decoupled from other fields prior to inflation. It has 
been shown in that this thermal distribution of inflaton modifies the power spectrum of 
inflaton fluctuations and the curvature power spectrum will have an additional temperature 
depended term. In this section we compute the two point correlation of inflaton perturba- 
tions taking this thermal distribution of inflatons into consideration. 

The Fourier expansion of inflaton fluctuations in de-Sitter space is 

/^3l, , . 
^(6k¥^,(t) + 6^„k¥'*,(^))e^''•^ (1) 
(27r)2 ^ ^ 

where (pk{t) are the mode functions which satisfy the Klein-Gordon equation in Fourier space 
and 6k and fej^ are the annihilation and creation operators respectively. In Fourier space the 
inflaton fluctuations can be written as 

6<l){k,t) = b^Mt) + bl^vm. (2) 

The canonical commutation relation satisfied by these creation and annihilation operators 
is 

= 6'{k^-k2), (3) 

with the vacuum satisfying 6k 1 0) = at zero temperature, which ensures that the vacuum 
has zero occupation Nk\0) = where = fe^^k is the number operator. The power 
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spectrum of inflaton P5,j>{k) is the two-point correlation function of the inflaton fluctuations 
in momentum space which is defined as 

Ps^(k)^^{6(l)ik,t)5<t>ik,t)), (4) 

where k = |k|. 

In this case, where there is a radiation era prior to inflation, the inflaton will have a 
thermal distribution during inflation. Due to this distribution the thermal vacuum \Q) = 
Irifci, nk2, ■ ■ ■) will now contain real particles yielding 

N,\Q)=nk\Q), (5) 

where is the number of particles with momentum k present in the thermal vacuum. In 
general, for creation-annihilation operators with different momenta one gets 

bib^M = SHk^~k2)nM. (6) 

Throughout this paper we will consider non-interacting real scalar field for which the chem- 
ical potential fi = 0. For a single inflaton with momentum k the partition function will 
be 

, = J2 e-^-^' = , (7) 

nfe=0 

where /3 is the inverse of the comoving temperature T. Due to this thermal distribution of 
the inflaton fluctuation a thermal statistical average of the two-point correlation function 
will determine the power spectrum 

1,3 

27r2 

^ J2pi^k){^mk,t)6<p{k,t)\n). (8) 



_ k^ 

~ 2^ 



Here p{ek) is the probability of the system to be in the state Ek = rikk which is defined as 



where z is given in Eq. ([7]). However, due to the thermal distribution of the inflaton field 
the inflaton fluctuations will follow the relations given in Eq. and Eq. ([5]) which yield 

{n\s<j){k,t)S(j){k,t)\n) = \^k{t)f{n\{i + 2Nk)\n) 

= \ipkit)f {1 + 2nk) . (10) 



Hence the power spectrum given in Eq. ([8]) will be 

fc3 



27r2 "^"^ " z 

Mt)\H^ + '2fB{k)), (11) 



2ti 
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where /b(A;) = ^jjfty Bose-Einstein distribution. To get the last equality in the above 

equation the following relation is used j27| 



oo 



E^^^ = 77^- (12) 

Now for a light infiaton (m^ -C -ff, being the mass of the inflaton and H being the 
Hubble parameter during inflation) the mode function has the solution 28 1 



l^^l - ^ f47 ) > (13) 



/2F \aH 

where a is the cosmic scale factor and z/^ ~ | — In a generic single field inflationary 
model this mode function solution along with the factor in the power spectrum gives a 
nearly scale invariant spectra for inflaton fluctuations. But due to the thermal distribution of 
the inflaton fluctuations, expression for power spectrum in Eq. flTTl) contains an additional 
temperature dependent factor of (1 + 2/b(A;)) = coth(/3A;/2). Thus the thermal power 
spectrum of inflaton fluctuations is given by 

PS(fc)=coth(/3fc/2)A<^(fc), (14) 

and hence the thermal average of the power spectrum for comoving curvature perturbations 
defined in Eq. flASP will be 

V'^{k) = coth{(3k/2)Vn{k), (15) 



as has been already stated in 17l |. 

In 17| the CMB power spectrum generated using the thermal comoving curvature power 
spectrum is compared with WMAP data and a constraint on comoving temperature is given 
as T < 1.0 X 10"'^ Mpc~^ (with the convention that the present scale factor = 1). 
Such a bound is also found in 1^ from thermal primordial gravitational waves. Since 
T = ajTph where Tph and are the physical temperature and the scale factor when our 



current horizon scale crossed the de-Sitter horizon during inflation, this constraint can be 
rewritten as Tq < A.2H. As the comoving wavenumber k = aiH one can put a lower bound 
on (3k as 

^k = ^— > 0.238. (16) 

'^i-' ph 

This lower bound ow (3k will be used in following sections to quantify the maximum value 
of non-Gaussianity in thermal bispectrum and thermal trispectrum. 

3. NON-GAUSSIANITY IN BISPECTRUM FROM THERMAL DISTRIBUTION 
OF INFLATON 

The three point correlation function of comoving curvature perturbations TZ or the bis- 
pectrum is deflned in Eq. f lAlOp and the non-linear parameter for bispectrum in the case of 



single fleld slow-roll model is given in Eq. ( 1A14I1 . In presence of a pre-inflationary radiation 



era the bispectrum will also receive a modification as in the case of the thermal power spec- 
trum. Hence in this case the three point correlation function of the non-linear curvature 
perturbation will be 

1 [hV d fH\ f d^p 



(7^A^L(kl)7^JVL(k2)7^^^L(k3))^ 2 



(2vr)^ 

(50L(p)50L(ki - p)50L(k2)50L(k3))^ + 2 perms] (,17) 

where R.H.S. of the above equation contains the thermal average of four-point correlation 
functions of the inflaton perturbations. 

We will flrst generalise the case of thermal average of the two-point correlation function 
to derive the thermal average of the four-point correlation function of scalar perturbations 
with any four momenta. The thermal average of higher order correlation functions is of the 
form 

(0fci0fc20fc3 ■ ■ ■ )/3 = Pih, h, h,---) (^I0fci0fc20fc3 ■ ■ ■ 1^), (18) 

where the thermal probability of the occupancy of different momenta kj and e = nk,.kr 
is 



r 



CO ^ 

HE = 1113^. (20) 



Here Z is the grand partition function of massless inflatons with energies E^^ = y = /c. 
which is given as 

oo 

z " 

where r is the index for different energy levels. 

The four-point correlation function of inflaton fluctuations with four different momenta 
contains six different combinations of two creation and two annihilation operators and ther- 
mal average of one of these combinations can be derived as follows : 

Let us consider the thermal average of ^&^ki^k2^-k3^k4 j which yields 

[bl^bu^bl^bk^^^ = ^pih,k2,h,h) (n^bl^^bk,bl^,bk^ 



5'{k^ + k2)5%ks + k4)^ 5^ 5^ ^-P{n.,k.+n,^k,) ^ ^21) 



Z 



where Z = (^jz^wi^ (r^^^^)' '^^^ summations in the above equation yields 

(b^A^bUA^)^ = ^'(ki + k2)5'(k3 + k4) [fBih)fBih)] , (22) 

where the identity stated in Eq. (fT2l) is used. Similarly the thermal average of other combi- 
nations of the two creation and two annihilation operators will be 

{bi.M.bUblu)^ = 5-'ik, + k,)6'ik, + k,)[l + Mh)] + 

S\ki + k,)6%k2 + k4) [1 + fBih) + fBih) + fB{h)fB{k2)] ,(23) 
(&k,6^„k,&k3&^_k,)^ = '^'(kl + k2)5^(k3 + k4) [1 + fBih) + fBih) + fBih)fBih)] ,(24) 

'Mbl^/_^^b^,) ^ = 6'ik, + k2)53(k3 + k4) [fBih) + fBih) fBih)] , (25) 



&Lk,&k,&k3&Lk,;^ = 5lkl + k2)(5^(k3 + k4) [fBih) + fBih) fBih)] , (26) 

(&U,^-k,&k3&k4)^ = -5'(ki + k4)5^(k2 + k3)/B(A;i) 

+5%k, + k,)5%k2 + k4) [fBih) fBih)] . (27) 

Hence the thermal average of a general four-point correlation function with four different 
momenta will be 

(50(ki, t)6<Pik2, t)50(k3, t)6<Pik,, t))^ = it)f WkM? [5'(ki + k4)5=^(k2 + k3) 

+5\k^ + k;)5\k2 + k4) {1 + fBih) + fBik2) + 2fBik^)fBik2)}] + |<^A,(t)|' WkM? 

X [5\k^ + k2)53(k3 + k4) {1 + 2/b(A;i) + 2/b(A;3) + 4/b(A;i)/b(A;3)}] • (28) 
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With this general result one can calculate the thermal average of the three-point correlation 
function of the comoving curvature perturbations using Eq. f|T7|) as 



X 



'Pnih 


) Pn{k2) 1 




kl \ 


Pn{k2) 


Pnih) ( 


1.3 


kl V 


Pn{h) 


Pniki) ( 


1.3 


kl V 



l + \fB{kr) + ]^fB{k^) + fB{k,)fB{k,] 



+ ^ ^ '\T' ( 1 + \fBik2) + lfB{k3) + fB{k2)fB{k3) 



(29) 



where Pn{k) is defined in Eq. (1A11I) . The three momenta form a triangle due to the presence 
of the delta function. In general three different triangle configurations are considered to 
determine the non-Gaussian effect. The non-linear parameter /tvl for these three momenta 
configurations are discussed below: 



Squeezed triangle case : For a "squeezed" triangle the configuration suggests |ki 
|k2| ~ A; ^ Ikgj. In this configuration the /at^ will be 



Inl = - e) (^2 + 2Mk,) coth i^^jj . (30) 

At low temperature /? — > 00 and f'Bik^) 0, yielding the same contribution to the f^L 
for super-cool inflation. The minimum value k^ can obtain when the corresponding 
wavelength is of Hubble size while crossing the horizon such that -^3 = ~ which 
implies pk^ ~ 0.238. Hence it yields 



Pnl = - e) X 2 (^1 + 3.72 coth 

= fNLX 2(^1 + 3.72 coth (^^^ ^ . (31) 

A lower bound on f3k can be given from thermal power spectrum which is given in 
Eq. ( fT6ll and for this constraint f^L will be maximum yielding f^j^ = 64.82/^1 ~ 0.65. 

Equilateral triangle case : For a "equilateral" triangle we have |ki| = |k2| = jkgl = k 
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and in this case the /nl will be 



fNL = g(^-e) (3 + 
= /^L I 3 



5 



4sinh2 (f ) 
5 



(32) 



4sinh2(f) 

This implies that for the modes corresponding to our present horizon f3k > 0.238 and 
the fj^L = ~ 0.9. 

Folded triangle case : For "flattened" isosceles triangle or the "folded" triangle case 
we have |ki| = [ksl = ||k2| = k and in this case the /nl will be 



fNL = liS-e)\3 + 



sinh^ (f ) 



NL 



sinh^ if) 



In this configuration the non-linearity will be / 



th 
NL 



(33) 



73. 28 f^L ~ 0.73 at horizon 



crossing for the modes corresponding to our current horizon. 
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FIG. 1: as a function of dk for different triangle confiffurations of the three momenta. 

JNL ' ° ° 

fth 

In Fig. ([1]) the thermal enhancement factor is plotted as a function of f3k for three 
different triangle configurations of the three momenta. From the above discussion it is seen 
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that the maximum contribution for Jnl comes from the "equilateral" configuration, though 
the contribution from the other two configurations are of the same order. Non-Gaussianity 
in all these three cases may be measurable by the 21-cm background radiation observations 

3- 



4. NON-GAUSSIANITY IN TRISPECTRUM DUE TO THERMAL DISTRIBU- 
TION OF INFLATON 

In a generic slow-roll single field model the trispectrum T{ki, k2, k^, k^) is defined as the 
Fourier counterpart of the connected part of four point correlation function of comoving 



curvature perturbation 



30| 



(7^(kl)7^(k2)7^(k3)7^(k4)), = TNLS^k, + ka + kg + k4)T(A;i, /c2, k^, k^), (34) 

where tnl is the non-linear parameter for trispectrum and is a measure of non-Gaussianity. 
The connected part in the above equation is defined as 

(7^(kl)7^(k2)7^(k3)7^(k4)), = (7^(kl)7^(k2)7^(k3)7^(k4)) 

- ((7^i(kl)7^L(k2))(7^i(k3)7^L(k4)) + 2perm). (35) 

The comoving curvature perturbation TZ has been expanded non-linearly upto 0(501). 
Hence the term (7?.L(ki)7?.L(k2)7^L(k3)7?.ArL(k4)) vanishes as it turns out to be a expec- 
tation value of odd number of Gaussian variables 50^. If one expands TZ upto (9((50|) then 
a term like (7^i(ki)7^2.(k2)7^L(k3)7?./vL(k4)) will survive and the non-linear parameter Tjy]^ 
will be ttvl = (I/atl)^ i-e. ((9(e^)) [14]. The trispectrum T{ki, k2, k^, k^} in this case will be 
proportional to product of three power spectrums and the four momenta form quadrilateral 
configuration due to the delta function in Eq. flM|) . 

But in this model of slow-roll inflation with a radiation era prior to inflation, the analysis 
for trispectrum turns out to be quite different. In presence of a pre- inflationary radiation era 
the four point correlation function which contributes to the non-Gaussianity will be thermal 
averaged as in the case of power spectrum and bispectrum. It is worth to point out that due 
to thermal averaging the four point function is not just the square of the two-point function 
as that would have been the case at zero temperature. So in this case, by connected part 
of the four-point function defined in Eq. (I35l) we will simply mean the excess of the thermal 
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averaged of four-point function than the square of its two-point Gaussian part and will now 
define the non-linear parameter t^l in the following way 

(7^i(kl)7^i(k2)7^i(k3)7^i(k4)), = (7^i(kl)7^i(k2)7^i(k3)7^i(k4))^ 

- ((7^L(kl)7^L(k2))^ (7^L(k3)7^L(k4))^ + 2 perm.) 

= -NL i^^5^(kl + k3)5^(k2 + k4) 

+2 perm.] . (36) 

Hence in this case tnl will not depend upon the slow-roll parameters. The thermal average 
of the four-point correlation function of inflaton fluctuation has been calculated in the last 
section in Eq. (1281) . Using this equation the thermal average of the four-point correlation of 
curvature perturbation can be derived as 

(7e4ki)7e4k2)7e4k3)7e4k4))^ = ^^l?'^ ^'^[^'^ [^'(ki + + k3) 



+5\k, + k3)<5^(k2 + k4) {1 + fnih) + fB{h) + 2fB{h)fB{k2)}] + ^""^ ^""A^'^ 
X [^^(k^ + k2)(53(k3 + k4) {1 + 2fB{ki) + 2fB{h) + ^fB{kl)fB{h)}] , (37) 

and the thermal average of two-point function can be given in terms of the power spectrum 
as 

(7ei(ki)7eak2))^ = ^^^^^ (1 + 2fB{k{)) 5\\i, + k2). (38) 
Hence the connected part will be 

(7^akl)7^z.(k2)7^ak3)7^z.(k4)), = -^-^^^^ [S%k, + ks)6%k2 + k,) {fB{k,) + 

^2 

fB{k2) +2fB{kl)fB{k2)} + 26^{k, + k^)6^k2 + k,) {fB{k,) + fB{k,) + fB{ki)fB{k2)}] . 

(39) 

The four momenta in this case will not form a quadrilateral as in other trispectrum cases. 
But due to the presence of two delta functions on the R.H.S. of the above equation the 
non-linear parameter t^l can be calculated in the following two cases : 

1. ki = — ks, k2 = — k4 and ki = k {i = 1, 2, 3, 4) : 

^^'^ = ~cosh(/5A;)-l- 
12 



The maximum observable value of \t^l\ '^^^ t)e obtained using the constraint on 
the comoving temperature as j3k > 0.238. Hence for f]k ~ 0.238 one finds that 



r*,^ ~ -35.14. 



2. ki = — k4, k2 = — ks and ki = k {i = 1, 2, 3, 4) : 

1 - 2ef"' 



th 

'NL 



-2- 



The maximum value of r^i^ for this case will be t] 



th 



NL 



-42.58. 



(41) 




FIG. 2: Plot of for two different momenta configurations as a function of Pk 



In Fig. ([2j) we have plotted as a function of pk. We find that the maximum value of 
non-Gaussianity comes from the configuration when ki = — k4, k2 = — ks and ki = k {i = 
1,2,3,4) which is ~ —42. We do not compare this contribution to non-Gaussianity due to 
thermal initial states with the zero temperate case as there is no contribution from the later 
at this order and the leading order tnl in zero temperature is 0{e^). 



5. CONCLUSION 



We studied the effect of a decoupled thermal spectrum of infiatons (which exist in the 
scenario where the inflation is preceded by a prior thermal era) on the non-Gaussianity of the 
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primordial perturbation. We found that thermal inflatons can enhance the bispectrum non- 
Gaussianity parameter /atl by a factors of (65 — 90) depending upon the momentum configu- 
ration. The zero temperature non-Gaussianity parameter /tvl in single field inflation models 
is proportional to the slow roll parameters and is expected to be of order ~ 10^^. Therefore 
the observed value of f^L in thermal history models will be of ~ 1. This is too small to 
be measured by WMAP or even the forthcoming PLANCK experiment. Measurements of 
anisotropics in the Hydrogen 21-cm radiation background can detect non-Gaussianities as 
low as /tvl ~ 0.1 9|, and this may be the ideal experiment in which non-Gaussianities with 
a thermal origin can be observed. The 21-cm observations may also be able to measure the 
non-Gaussianity in the trispectrum Ttvl ~ (9(10) and for which the prediction from thermal 
history inflation scenarios is > t/vl > —43. We conclude that a signature of thermal infia- 
ton background at the time of inflation is a large trispectrum non-Gaussianity compared to 
the bispectrum non-Gaussianity. 



APPENDIX A: NON-GAUSSIANITY IN A SINGLE-FIELD SLOW-ROLL IN- 
FLATION MODEL 

In a single field slow-roll model Non-Gaussianity appears generically once the inflaton 
field has self-interactions such as V{(f)) ~ 0^ or V{(f)) ~ 0^. But the Non-Gaussianity 
due to these self-interactions is very small where the bispectrum Non-Gaussian parameter 
/nl ~ C^(£^)- Larger contribution to Non-Gaussianity in such models come from non-linear 
curvature perturbations and in this cases /nl ~ C>{e,ri) [29]. Here we will briefly discuss 
how non-linearity in comoving curvature perturbation gives rise to Non-Gaussianity of the 
order of slow- roll parameters e and r]. 

The quantum fluctuations in inflaton field generates fluctuation in the metric which is 
coupled to it through Einstein's equation. The perturbed FRW metric has the form (con- 
sidering only scalar perturbations) 

[l + 2A A 

V -(1 - I 

where the quantity ip is known as the spatial curvature perturbation, r] is the conformal time 
and a{ri) is the cosmic scale factor. The gauge invariant quantity formed out of this spatial 
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curvature perturbation is known as the comoving curvature perturbation and defined as 

n{t, x) = tpit, x) + -^6(f){t, x). (A2) 



These perturbations are conserved on super-horizon scales throughout the evolution. Hence 
during inflation in the spatially flat gauge this reduces to 

7^(^,x) = :|50(^,x), (A3) 

and after inflation when 50 ~ this represents the gravitational potential on comoving 
hypersurfaces 

7^(^,x) = V^(t,x). (A4) 

The CMB anisotropy spectrum is determined by the power spectrum of this comoving cur- 
vature perturbation which is related to the power spectrum of scalar perturbation as 

P 1 

T'nik) = —{n{k)n{k)) = ^/sAk), (A5) 

where the slow-roll parameter e = AnGj^z, mpi = "^7^5 is the reduced Planck mass and G 
being the Newton's Gravitational constant. 

Presuming that the inflaton fluctuations S(f) are initially Gaussian, the comoving curvature 
perturbations TZ given in Eq. flASP also obeys Gaussian statistics in the linear order 

7^i(^,x) = ^(50L(t,x), (A6) 



where 7^L()f:,x) and can be expanded in Fourier space as 

7^i(^,x) = / ^e^^■''7^i(^,k). (A7) 

J (27r)2 

Being constant in time outside the horizon these comoving curvature fluctuations after en- 
tering the horizon in later times produces curvature perturbations which are Gaussian in 
nature. 



In the non-hnear limit one observes that 29] -j = — ^y^u) ^ function of and hence 



ml^ V"(0) 

7e^i(t,x) = S-6Mt,^) + (j) ¥i(i,x) + 0{6<j>l). (A8) 



15 



Therefore in the Fourier space one gets 

UNLit, k) = -j6Mt, k) + T^TTT T / -^^Mt, p)SMt, k - p), (A9) 

where S(f)L{t, k) has the same form as given in Eq. ([2]). 

The non-hnear parameter f^^ for bispectrum B{ki, ^2, k^) or the three-point correlation 
function of the comoving curvature perturbation is defined as jsoj 



(7^(kl)7^(k2)7^(k3)) = (27r)-^(5^(ki + k2 + k3)-/jVL + 2 perms. 

= (27r)-i53(ki + k2 + k3)^/^ii?(A:i, A^s, fcs), (AlO) 

5 

where we have used the definition 

Pn{k) = {2n^)Vnik), (All) 

and the normahzation (27r)~t has been chosen accordingly. Hence using Eq. (\A9h one can 
compute the bispectrum in this follows 

(7^^i(k07^^L(k2)7^;vL(k3)) ~ i2n)-l5%k^ + ks + k3)(2m^ie)^ 



V~fcf kf~ J ■ ^ ^ 

Comparing the above two equations the non-linearity parameters in this case will be 

5 ^ d fH 

~ 3 90U'(0)J- ^ ^ 

The non-linear parameter /tvl given in Eq. flAlSp can be fully expressed in terms of the 
slow-roll parameters e = (v')^; V = (t^) ^^'^ 6 = rj — e as 

fNL = ks-e). (A14) 



The form of the f^i derived above is same as in [1^]. Here f^i being proportional to the 
slow-roll parameters is too small to be detected by the ongoing experiments. 

The non-linearity parameter can also be expressed in terms of the potential V{(f)) using 
Eq. flAT3l) as 



5 , (V'\ d fVi<P)\ 
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This equation is useful in the case where the form of the potential is known. Such as, for a 
power-law potential where V{(j)) ~ 0" the non-linearity parameter will be 



NL 



5 \ m 
- n 
6 



2 

PI 



(A16) 



which is same as the one given in 



29|. 
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